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Abstract 

We study a kind of filtering, an amplitude truncation with upper and lower truncation levels 
i max and x m i n . This is a generalization of the simple transformation y(t) = sgn[x(£)], for which a 
■ rigorous result was obtained recently. So far numerical experiments have shown that a power law 

U: spectrum 1//" seems to be transformed again into a power law spectrum 1/ f@ under rather general 

condition for the truncation levels. We examine the above numerical results analytically. When 
1 < a < 2 and x max = —x m i n = a, the transformed spectrum is shown to be characterized by a 
certain corner frequency f c which divides the spectrum into two parts with different exponents. We 
derive f c depending on a as f c ~ <x _2 '( a_1 ). It turns out that the output signal should deviate from 
the power law spectrum when the truncation is asymmetrical. We present a numerical example such 
that l// 2 noise converges to 1// noise by applying the transformation y(t) — sgn[x(t)] repeatedly. 
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1// noise whose power spectrum is inversely proportional to the frequency / has been observed in a 
variety of systems, since it was first discovered in current fluctuations of a vacuum tube (see ref. 1 
and references therein). Sometimes l/f a noise with the exponent a not close to 1 but between and 
2 is also referred to as '1// noise'. When this extended definition is applied, '1// noise' is more 
widely observed.!) 

It should be noted that the integral f£° f~ a df diverges for all a. Therefore, if the process is 
stationary, which is usually a reasonable assumption, a low or high frequency cut off should exist 
corresponding to a > 1 or < a < 1, respectively. Or equivalently, the value of a cannot be constant: 
a should become smaller (larger) than unity in a low (high) frequency range. The power law behavior 
involves several difficulties. The autocorrelation function of 1/ f a noise with a < 1 decays very slowly 
(power law decay, see eq. (2.10b) and Fig. [j](a)) especially when a is close to unity, unless there is 
another low frequency cut off /o so that the spectrum becomes white (a = 0) in a range / <C /q. 
It often happens that the cut off frequency /n is extremely low or not observed during the practical 
observation timeE'B One of the difficulties is to explain a long time scale ~ l//o from any realistic 
model. Another problem is the value of the exponent which largely deviates from 2. The exponent 
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a = 2 is most naturally expected in nonstationary processes. Random walk in n-dimensional space 
yields 1/ f 2 spectrum. Also, a Lorentzian spectrum (/ 2 + /d) _1 looks as l// 2 when the observation 
time is not long enough: / S> /o- Lorentzian spectrums are obtained from Debye-type relaxation 
processes, which occur most commonly. A randomly amplified Langevin system which has a power law 
distribution function also leads to a Lorentzian spectrum^ even though the process is nonstationary 
so that the Wiener-Khintchine theorem does not hold. Then what mechanism could exist that would 
reduce the value of a from 2? The present work relates to this second problem. We shall investigate, 
by generalizing a previous theory ,0* how Gaussian 1/ f a noise is affected by a kind of filtering, an 
amplitude truncation. 

In the previous paper, a simple dichotomous transformation defined by y(t) = sgn[x(i)] for Gaussian 
l/f a noise was studied in detail and a rigorous result for transformation properties was obtained.0* 
When a Gaussian noise x(t) with a power law spectrum 1/ f a is filtered by this transformation, power 
spectral density (PSD) of the output noise y(t) obeys again a power law The exponent [3 is 

derived as (3 = a for < a < 1, and (3 = (a + l)/2 for 1 < a < 2. 

The present paper deals with a more general amplitude truncation 

!^min if X\t) ^ 3?min> 
£max if x(t) > X max , (1.1) 
x(t) otherwise, 

for Gaussian l// a noise. The transformation y{t) — sgn[x(£)] is equivalent to the situation x max — 
—Xmin <C 1 in eq. ( |1.1[ ), namely, the dichotomous transformation is a special case of eq. (|1.1[). So far 
numerical experiments have showrJli) that a power law spectrum 1/ f a seems to be transformed again 
into a power law spectrum l/f@ under rather general condition for the truncation levels x m i n and 
a^max- However, the numerical results are not always clear, so that we will investigate analytically the 
above amplitude truncation to know how (3 depends on the levels as well as a and when the output 
signal deviates from the power law spectrum. 

The rest of the present paper is organized as follows. In the next section, we briefly review and 
reinterpret the recent result for the dichotomous transformation. In §3, we extend our study for the 
symmetrical truncation with finite levels : x max = — x m i n = a in (|1 . 1|) . We consider the asymmetrical 
truncation in §4. It will be shown how PSD of the output signal deviates from the power law. In §5, we 
present a numerical example such that 1 // 2 noise converges to 1 // noise by applying the dichotomous 
transformation repeatedly. The last section is devoted to summary and discussions. 



§2 Dichotomous Transformation 

The dichotomous transformation is defined by 

9 ( ( , =sg n W )] = {;; I 'BSJ; (2.1) 

where x(t) is an input noise with zero mean and y{t) is an output noise. This is a special case of 
the amplitude truncation (1.1) and was investigated in detail in ref. 1. It turned out that when the 
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transformation fl2.1| ) is applied to a Gaussian l/f a noise x(t), then the output noise y{t) exhibits a 
power law spectrum 1/ f@. The exponent f3 was derived as (3 = a for < a < 1, and /3 = (a + l)/2 
for 1 < a < 2. As will be shown below, a key to understand this result is a transformation property 
of the correlation function R(t) (eq. Q2.7|) ). 

The correlation function of the output signal y{t) is obtained from eq. (|2.1|) as 



R y {t) EE (y(t )y(t + t)) 

= IP (y(O)y(t) = 1) + (-1) • P (y(O)y(t) = -1) 
= P (x(0)x(t) > 0) - P (x(0)x(t) < 0) 

= 2 ■ P (z(0)a:(t) > 0) - 1 , (2.2) 

where P(.) stands for the probability that the condition of the argument is satisfied. In the above, 
we have assumed a stationary process for the input noise, which leads to a stationary process for 
the output one. We can therefore apply the Wiener-Khintchine theorem to both input and output 
noises. We also assume that the process for the input noise is correlated Gaussian.!) That is, the joint 
probability is given by 



P(x(t ),x(t + t)) = P(x(0),x(t)) 
= -rexp[-(x 2 
ee f(x,y,c) , (2.3) 



^-exp[-(x 2 -2cxy + y 2 )/B] 



where, 

x ee x(0) , 

y = x{t) , 

c ee R x {t) , (2.4) 

A ee 2irR x (p)y/T=cl , 
{ P = 2P,(0)(l-c 2 ) , 

and 

Rx(t) = (x(to)x(to+t)) (2.5) 

is the correlation function of x(t). 

The probability P(xy > 0) is calculated from eq. (p.3[) as 



P{xy>0) = dx dyf(x,y,c) 
Jo Jo 

+ dx dyf(x,y,c) 

J—oo J —oo 



= - + - axcsin^t)] . (2.6) 

Thus we obtain the relation between the input and output correlation functions as 

2 

Ry(t) = -arcsinfP^t)] . (2.7) 

7T 
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The previous result for the exponent (3 can easily be understood qualitatively if we note the approx- 
imated expression for (|2.7|): 

when \R x (t)\ < 1 , R y (t) oc R x (t); (2.8) 
when |1 - R x (t)\ < 1 , Ry(t) oc const. - (1 - R x (t))? . (2.9) 




Fig. 1. Correlation function R x (t) of a l// a noise, (a) < a < 1: eqs. ( 2.10 ). (b) 1 < a < 2: eqs. 

O- 

First we explain /? = a for < a < 1. When PSD for low frequencies obeys the power law as 
~ f~ a with < a < 1, the corresponding correlation function R x {t) should be R x (t) ~ for large 
t values.00) Let us consider a system whose characteristic time scale is t\. Then, by assuming R x {t) 
as 

R x (t) = 1 for t/n < 1 , (2.10a) 

= (*M) a_1 for t/n > 1 , (2.10b) 
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the PSD is ~ (ti/) _q for low frequencies: < / <C 1/tl. The relation ( |2.8| ) can be used to obtain 
Ry(t) for large t: t/r\ 3> 1, because R x (t) is small. This is why the same exponent (3 = a is obtained 
for < a < 1. 

The case 1 < a < 2 can also be understood qualitatively. Let us assume the correlation function as 

R x (t) = \-{t/Ti) a ~ l for t/n<l, (2.11a) 
R x {t) = for t/n > 1 . (2.11b) 

The above correlation function yields a white spectrum if the observation time t is long enough: t 3> tl, 
i.e., / <C 1/ti. Naturally we are interested in the situation in which the characteristic time scale t\ 
is very long. The power law, eq. ( p. 11a ) again leads to the power law spectrum S x (f) oc (ti/)~ q for 



a wide range: / 3> 1/nxP Note that if t\ is very long, the lower limit 1/n can be small and the 
power law ~ / _Q holds practically in a wide range of low frequencies. When time t is long enough 
but satisfies t <C T\, the approximation Q2.9| ) can be used for R y (t): 

a — l 

R y (t) oc const. — (t/ri)^~ 

= const. - (t/ri)^ 1 . (2.12) 

Corresponding to the above R y (t), the PSD is obtained as ~ f-P, where (3 = (q + 1)/2. 

§3 Symmetrical Truncation with Finite Levels 

In this section we generalize the results obtained in the above. Suppose that a Gaussian l/f a noise 
with zero mean is transformed by a symmetrical truncation with finite levels ic max — — ^min = 

Let 

us define a typical time scale t c in which the noise signal passes between the two truncation levels. 
In a high frequency range / 3> f c = 1/t c , PSD of the truncated noise is mainly determined by the 
behavior of the noise between the two levels. On the other hand, in a low frequency range / < / c 
(corresponding to a long-term observation), PSD of the output signal should have a similar PSD to 
that obtained by the dichotomous transformation. Therefore, for < a < 1, the same exponent (3 = a 
is expected both in high (/ c < / <C 1/ti) and low (/ < / c ) frequency ranges. 

In contrast, when 1 < a < 2, PSD is divided into two parts by a corner frequency / c : = a for 
a high frequency range (/ > / c ) and (3 = (a + l)/2 for a low frequency range (1/ti <C / < / c )- A 
typical PSD of a truncated signal for a = 2 is shown in Fig. |2[ 

The aim in this section is to derive the dependence of f c upon the truncation level a. The transfor- 
mation is expressed as 

la- sgn[x(t)J otherwise, 
where x(t) is a correlated Gaussian noise with zero mean. 

As in the previous section, the joint probability (|2.3|) yields the correlation function of the output 
signal as 

R y (t) = 2 .a 2 (I(c)-I(-c)) + 2.(J(c)-J(-c)) 

+A-a{K(c)-K{-c)) , (3.2) 
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Fig. 2. PSD of the truncated signal for a = 2. The PSD is obtained by averaging over 40 samples. 
PSDs in three frequency ranges are composed. Straight lines are S{f) oc f~ 2 and oc / _1 ' 5 



where 



/(c) 



OO POO 



dx I dyf(x,y,c) , 



pa pa 

J(c) = / / dyxyf(x,y,c) , 
Jo Jo 

roc pa 

K(c) = / dx dyyf(x,y,c) . 

Ja JO 

If we assume oCl, /(c), J(c) and K(c) can be expanded with respect to a as 
/(c) = -\l + - arcsin(c)] - - a 2 ^ + 0(a 3 ) , 

4 7T t4q 



J(c) 
K(c) 
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+ C»(a 5 



2 vBtt o2c - 3 . 4 . 

+ a 3 - , + Q(a 4 ) . 



2\2A 3lA 

Substitution of Eq. (E01) into (|3. 2|) leads to the correlation function of the output signal as 



(3.3) 



R y (t) = 2a 2 {^ arcsin(c) - 2 • 



ca 



27TV1 - c 2 



+ •••} 



+4a{ 



ca 



+ •••} 



a 2 {— arcsin(c) 



' 7T 



7T 



2ca 2 4ca 2 

+ — -7=^ + -"} 
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Since we are considering the case of 1 < a < 2 , we substitute c = R x (t) = 1 — t a 1 (eq. fl2.11a| )) into 
the above. The correlation function is thus obtained as 

R y (t) ~ a 2 {l - ^?t(«-!)/ 2 



7T 

a-1 



-a 2 + •••}, (3.4) 

for t <C 1. Here and hereafter, the dimensionless time i/ri is replaced by t for simplicity. Thus t<Cl 
actually means i < ti. If a is small enough and/or t is not so small so that the third term in the 
braces in (3.4) can be neglected, R y (t) reduces to eq. ( |2,12j ). As t decreases, the third term grows 
and becomes comparable to the second term. In other words, the perturbative expansion (3^4) breaks 
down when f( a-1 )/ 2 ~ a 2 t~^ a ~ 1 ^ 2 i.e. t ~ a 2 '^ -1 ). This means that the corner frequency f c depends 
on the truncation level a as 

/c-a- 2 ^- 1 ). (3.5) 

For a long-term observation, i.e., 1 >t > a 2 /( Q_1 ) (which corresponds to 1 -C / < f c ), the exponent (3 
of the PSD becomes j3 = (a + l)/2 because ( |3.4|) reduces to ( 2.12j ). For a high frequency range: / > / c , 
the same exponent (3 = a as the input signal is expected by the reason mentioned already, although 
the correlation function of the form R y {t) ~ const. — i"" 1 (which corresponds to the PSD~ f~ a ) 
cannot be obtained from any correction terms added to Q3.4| ) because the perturbative expansion has 
broken down. One can see in Fig. Q that the PSD obtained from the numerical simulation is really 
composed of two power law spectra separated by the corner frequency f c . 



The scaling property of eq. (3.5) can also be suggested from the following self-affine character of 
the signal from the fractional Brownian motion: 



{x{tf) ~ t 2H . (3.6) 



It is well known that if a signal satisfies the scaling relation ( |3.6| ), then its PSD is 1/ f a with a = 2H+1 
where a > 0. The above self-affine character yields 

-yj(x(ahy). (3.7) 



oc 

a 



This means that if the truncation levels are changed from ±J7 to ±aU and the time is changed from 
t to a l l H t simultaneously, the structure of the noise is invariant. Then the corner frequency changes 
from f c to a~Q' H ' f c = a _2 '^ a_1 )/ c . This coincides with the result obtained from the perturbation 
method; the present argument can be applicable even for large a values. 

To confirm a-dependence of f c , eq. (|3.5| ), we performed a numerical simulation for a = 2 by varying 



the truncation level a in eq. (3.1). For each a value, a PSD S(f) was obtained by averaging over 40 
samples. Then log S(f) was fitted to the line F(f) = (-1.5 log / + 1.5 log f c + A) 6 (log f c - log /) + 
(-2.0 log / + 2.0 log / c + A)6(log/- log / c ) where &{x) is a Heaviside function with 9(0) = 1/2. The 
result is shown in Fig. |3|. Agreement between the theoretical line: f c oc a~ 2 and the numerical result 
is tolerable. 
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Fig. 3. Log- log plot of the corner frequency f c vs the truncation level a for a = 2. The straight line 
is f c oc a~ 2 . 

§4 Asymmetrical Truncation 

We consider the simplest case, i.e., the following asymmetrical dichotomous transformation: 

m=^w) -.]={:; <«) 

The input signal is a Gaussian l// a noise with zero mean, and the condition < a < 1 is assumed 
here. 

Figure |I| illustrates the PSD of the output signal obtained by eq. ( |4.1| ) for a = 0.8 and a = 0, <r, 2<r, 
where <r is the standard deviation of the input signal. It is unclear from these numerical results 
whether the PSD obeys the power law or not. 

Note that the mean value of the output signal is finite: 

/DO 
sgn[x - a]f(x)dx 
-oo 

= -2 r f(x)dx = y , (4.2) 
j o 

where f{x) is the Gaussian distribution function for the input signal: 

= ^2^(0) 

The correlation function of is thus given by 

R y (t) = ((y(to)-y)(y(t + t)-y) 



2R X (0) 



(4.3) 



S 



o 

1—1 

O 




Fig. 4. Output PSDs obtained from the asymmetric transformation ( [4,l[ ) where the exponent of the 
input noise is a = 0.8. Each PSD was obtained by averaging over 40 samples, (a), (b), and (c) 
correspond to the truncation level a = 0, a and 2a, respectively, where a is the standard deviation of 
the input signal. 



= (y(to)y(to + *)) ~{yf 

= 2P((x(0) - a)(x(t) - a) > 0) - 1 - (y) 2 . 
The probability P((x(Q) — a)(x(t) — a) > 0) is expressed as 

P((x(0) - a)(x(t) - a) > 0) 

/*oo ra poo roc -. 

2 1 dx / dy + 2 / / f(x,y,c) , 



where f(x,y,c) is the joint probability defined by (|2.3|) ~ (|2.5|) . The first term in 
but not on t: 



2 I dx dyf{x,y,c) =2 dxf{x)dx = —y , 

-co JO JO 



(4.4) 
(4.5) 



(4.6) 

depends on a 
(4.7) 



which yields 



P((x(0) - a)(x(t) - a) > 0) = -y + T , 
T = 2 dx dyf(x,y,c) . 

J a J a 



(4.8) 
(4.9) 



We calculate R y (t) in two cases a <C 1 and a 3> 1. Since i? x (t) is assumed as eqs. ( 2.10 ), we set 
Rx{0) = 1 hereafter. 
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When a <C 1, y is expanded with respect to a: 

2a 

y = ■ 



+ 0(a 3 



The second term T in (W^) is similarly expanded as 



m 1 1 

T = — I — arcsm c ■ , 

2 7T 11 TrVl^ 



2 1 - c 2 

a H = a +••• . 



(4.10) 



(4.11) 



For t > 1 the correlation function becomes c = R x (t) = t a < 1. Then T can be approximated as 



T=\- -La + -c + 1(1 - c+ ^c 2 )a 2 + 0(c 3 ) . 

2 V27T VT 7T 2 



Substituting (gl|, (gig ), and c = t^ 1 into fl£ID , we obtain 



Ry(t) 



2{-y + T)-l-y 2 



1(1 - a 2 )c+ — c 2 

2 



I 2 
7T 



-(l-a / )t Q! - i + — t 

7T 7T 



2(a-l) 



for i > 1 and a < 1. 



(4.12) 



(4.13) 



(4.14) 



Equation ( [4.14 ) implies that the correlation function of the output signal obeys a power law ~ t a ^ 
with a small correction term. When the first term in ( 4.14 ) is dominant, i.e. a 2 t a_1 -C 1, the PSD of 
the output signal turns f~ a . The condition a 2 t a ~ l <C 1 corresponds to / <C /i, where 

/! = (l/a 2 ) 1 ^ 1 -") . (4.15) 

The frequency f\ can be very large for small a values because 1 — a > 0. The power law spectrum 
/ _a is then observed in a wide range of frequencies. 



For a > 1, an asymptotic expansion of (|4.2[) yields 



-1 + 



cxp 



a 
~~2 



.{i + o ( i)} . 



(4.16) 



We can calculate T as follows. 



dxdy 



exp 



:r 2 — 2cxy + y 2 



1 



OO 1 







-exp 

2Y2 



1-c 2 
B 



Y 



a-cVY 



dX exp 



X 2 



Integrating by parts with respect to the variable Y, we obtain 
T 
E 



E+F+G , 
S 1 



(4.17) 



(4.18) 



G = 



A Q (1 - c 2 ) a 
B 

A)(i-c 2 ) 
5 



l-C 2 o 


POO 


X 2 ' 


~ B a _ 


/ dX exp 

J(l-c)a 


B 



-1 

aY r exp 

2Y2 



1-c? 



B 



-Y 



A (l-c 2 )J a 2 2Y 



dY —— exp 









exp 



_ dX exp 

a-cVY 

(a - cVY) 



X 2 



B 
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The second term F in ( 4.1S| ) can be neglected because a > 1. The asymptotic expansion of E is 
obtained by substituting X 2 = x and integrating by parts as 



E 



B 2 



2Aq(1 - c 2 )(l - c)a 2 



exp 



2a 2 (l - c) 



B 



O -j exp 



2a 2 {l-c) 
B 



(4.19) 



In a similar manner G is calculated, by substituting y = (\Y — ca) and integrating by parts, as 



G 



cB 2 



2A (1 - c 2 )(l - c)a 2 
Equations ( gig ) and ( |O0j ) yield 

m (!+ C ) 2 



exp 



2a 2 (l-c) 
B 



+ lexp 



2a 2 (1 - c) 
B 



(4.20) 



7ra- 



: exp 



1 + c 



+ o 



-^exp 



1 + c 



(4.21) 



Substituting JO|) , (gig ), O) and c = i"" 1 into flU), we obtain 



7ra z 



(1 + c) 2 



exp 



1 + c 



exp [—a 2 



(for a 2 > 1) 



(1 + 2t«-i + |t 2 («-D) eX p[-a 2 (l - T" 1 + t 2 ^" 1 )] - exp[-a 2 ] 



(for c«l). 



When a 2 t a_1 < 1, eq. ( |4.22| ) reduces to 



i? y (t) t^-^xpl-a 2 ] 



7T 



(4.22) 



(4.23) 



The above result leads to a 1/ f a spectrum for the output signal in a low frequency range: / <C /i, 
where /i is defined by ( 4.15[ ). In contrast to the case a< 1, however, f\ is small for a> 1. This means 
that the power law spectrum l/f a cannot be observed unless the observation time t is large enough 
to satisfy a 2 t a ~ l <C 1. When this condition is not satisfied, R y (t) includes correction terms which are 
£2(a-i) _|_. . Therefore PSD of the output signal deviates from the power law as A/ f a + B / f 2 ^ 1 + ■ ■ ■. 
Since a > 2a — 1, the second term makes the exponent of the PSD decrease. As the value of a or 
/ increases, the higher order terms become more important, namely, the exponent tends to decrease. 
Figure. || indicates this trend: The slope of the PSD turns flatter as a or / increases. 

§5 Construction of 1/ f Noise from 1/ f 2 Noise 

We present a numerical example to generate a 1// noise from l// 2 noises using the dichotomous 
transformation introduced in §2. The procedure is as follows: 1) Prepare an ensemble of independent 
Gaussian 1/ f a noises. 2) Transform these noises to dichotomous ones via eq. ([2.1|). As described in 
§2, PSDs of the transformed noises are proportional to 1//( Q+1 )/ 2 . 3) Average over a bunch of these 
dichotomous noises to re-Gaussianize, and prepare a number of such Gaussian l//( a+1 )/ 2 noises. 4) 
Go back to 2). 
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We applied the above method to Lorentzian noises with a long relaxation time t\, whose spectrum 
is ~ [f 2 + (1/ti) 2 ] -1 . The spectrum looks like 1/ f 2 when the observation time t is not long enough: 
t <C t\ so that / 3> 1 /t\ . The Lorentzian noises were generated by the following first order Markovian 
process 

1 " 

n. 



x(t + 1) = px(t) + n(t) ; p = exp 



(5.1) 



where n(t) is white noise with zero mean. As p approaches unity, the time series tends to a l// 2 noise. 
A typical spectrum obtained from ( |5,1| ) is shown in Fig. ||(a), where p = 1 — 0.001, i.e. t\ ~ 1000 and 
n(t) is a uniform random number in [—0.5,0.5]. Starting with such Lorentzian noises, we obtained 
the first and the second transformations as Fig. |B|(b) and Fig. ||(c), which are very close to 1/f and 
l// 1,25 as was expected. It is rather amazing that only five independent output signals were averaged 
to obtain each approximately Gaussian noise. 





1 r 









-1 


Q 


-2 


So 


-3 


o 

be 


-4 - 


o 






-5 




-6 




-7 





(a) 


i %\ 




\ 



















(b) 
















T 















(c) 













































-5 -4 -3 -2 -1 -5 -4 -3 -2 -1 -5 -4 -3 -2 -1 
logic/* logic/ \0g 1Q f 



Fig. 5. PSD of the signal in each step in the procedure to achieve a 1/f noise. PSD is obtained by 
averaging over 10 samples. The straight lines represent f~ 2 in (a), /~ L5 in (b) and / -1 ' 25 in (c). (a) 
PSD of a Lorentzian noise, (b) PSD of the noise transformed once, (c) PSD of the noise transformed 
twice. To obtain a Gaussian input noise, five output noises by the previous transformation were 
averaged. 



There are several methods to generate l/f a noise&E^lialiS ; among them, McWhoter's model 
for 1/f noise is well known because the model corresponds to real physical systems, for example, 
noise due to surface traps of carriers in a semiconductor© McWhoter's theory tells that if a large 
number of Lorentzian noises with various relaxation time r are superposed with an appropriate weight: 
w(t) oc 1/r in a range T2 < r < t±, the spectrum becomes 1// in the frequency range 1/t\ <C / <C l/Tg. 
In contrast to McWhoter's theory, the present method requires only one kinds of Lorentzian noises 
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specified by a single relaxation time t\\ w(t) = 5(t — n). Then the dichotomous transformation, if 
applied repeatedly, can lead to nearly 1// noise in the range / 3> 1/tl . 

In many cases, model systems exhibit 1// noise only when the system parameters are tuned to some 
special values. For example, simulations of phonon number fluctuations based on a Fermi-Pasta-Ulam 
lattice resulted in a 1// spectrum when the lattice size is N = 8.00 For larger N, however, the 
spectrum tended to Lorentzian ~ (/ 2 + /o) , although the time scale t\ = l//o became larger as N 
increases© In contrast, tuning is unnecessary in the present method of generating 1// noise, which 
is an essential point of the method. 

§6 Summary and Discussions 

We have analyzed the amplitude truncation of l/f a noises (eq. (1-1)). Although l/f a noises with a 
between and 2 are commonly called '1// noise', the transformation property is different depending 
on whether a is larger than unity or not. 

For a > 1, the output PSD under the symmetrical truncation obeys the power law with a smaller 
exponent (3 = (a + l)/2 (which is still larger than unity) in a low frequency range (/i -C)/ <C f c - 
The corner frequency f c depends on a and the truncation level a as f c oc a~ 2 ^ a ~ l \ Let us assume 
that the time scale of the system, n = is long enough so that the interval [fx, f c ] is sufficiently 
wide. Then the filtering effect in the measurement makes it possible to observe a power law spectrum 
with the exponent much smaller than 2, even if the original signal is close to Lorentzian: a ~ 2. As 
was shown by the numerical simulation, if we start from an ensemble of Gaussian l// 2 or Lorentzian 
signals and apply the symmetrical truncation and re-Gaussianization procedure repeatedly, the output 
signal converges to 1// noise. 

When a < 1, in contrast, the symmetrical truncation leads to a power law spectrum with the same 
exponent a. (A proof for a = 1 has been given only for the dichotomous transformation!!])) This 
means that the exponent less than unity can never be reached by the symmetrical truncation as far as 
we start from the exponent larger than unity. On the other hand, when the asymmetrical amplitude 
truncation is applied to l/f a signal with a less than unity, the output PSD deviates from the power 
law. The correction terms make the slope of the spectrum flatter as / increases. However, as far as the 
results of the numerical experiments (a = 0.8 in the present paper and a = 1 in ref. 5) are observed, 
the output PSD looks approximately as 1 / /"-like with an exponent smaller than that of the input 
signal. 

Amplitude truncation may generally occurs in systems composed of threshold elements, for example 
neural networks. It also occurs in the flow of packets in Internet systems where the overflow of 
the packets is deleted. Furthermore, measurement or analysis of signals often involves amplitude 
truncation. For example, a dichotomous transformation is used to analyze the (spatial) long-range 
correlation in DNA© If a Lorentzian or 1/ f noise is superposed on the original signal, filtering with 
such a truncation may easily leads to the observation of a l/f a noise. 

Throughout the present paper, we have assumed that the time scale of the system is long enough so 
that the original data exhibits the power law spectrum in a wide range of frequencies. As mentioned 
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already, to derive an extremely long time scale from any realistic model is a difficult problem. Once 
a long time scale has been assumed, however, the present study suggests a possible mechanism to 
generate l/f a noises with the exponent smaller than that of the original signal. 
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